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ABSTRACT 


An adaptive coherent multi frequency radar system has been developed for several app- 
lications. By means of this, the velocity distribution (Doppler spectrum) and spec- 
tral intensity of 15 different irregularity scales (waves and turbulence) can be 
measured simultaneously. Changing the azimuth angle of the antennas at regular inter- 
vals, the directivity of the wave/turbulence pattern on the sea surface can also be 
studied. 

Using this radar system, series of measurements for different air/sea conditions have 
been carried out from a coast-based platform in Southern Norway. Experiments in the 
Atlantic were also performed with the same equipment making use of the NASA Electra 
aircraft. 

There are many air/sea phenomena which play a role in relation to backscattering of 
radio waves in the microwave region: gravity waves "modulate" the capillary wave 

structure, overturning wave crests produce focussing effects and also periodic regions 
of strong turbulence, the boundary layer wind field with strong turbulence amplified 
by the ocean waves will conceivably leave a patchy and even periodic footprint on the 
sea surface. 

By virtue of the fact that our multifrequency radar allows us to measure the velocity 
distribution ("coherent and incoherent component") associated with 15 different ocean 
irregularity scales simultaneously in a directional manner, it is possible to study 
the different air/sea mechanisms in some degree of detail. 


1 INTRODUCTION 

Radio methods have a substantial potential in the study of the irregularity structure 

of the sea surface for several reasons: 

- They provide a rather unique possibility to measure the directional wave spectra 
with good directional resolution 

- By means of a raultifrequency radio method wavelength and wave velocity can be mea- 
sured independently 

- It is also possible to distinguish between coherent wave motion and incoherent 
motion ( turbulence ) 

We shall present a detailed theoretical and experimental study of the sea surface. 

We shall in particular demonstrate that a microwave illuminator can be tailored so 

as to optimize the coupling between electromagnetic waves and ocean waves. 

There are three different problem areas: 

A. One must obtain an understanding in regard to the coupling mechanisms between the 
electromagnetic waves on the one hand, and the sea surface on the other. 


179 



B. 





Having established information about the sea surface as a scattering surface for 
radio waves (the delay function f(r,t) or the spectrum E(K,<jd)), it remains to 
establish the relationship between this "description domain of the radio scien- 
tist" and that of the ocean scientist, who want information about most signifi- 
cant wave-height, wave-height spectra etc. 

To convert from scattering cross-section to wave-height one needs information 
about how gravity waves and other large scale phenomena affect the small scale 
phenomena (capillary waves) which are comparable with the wavelength of the radio 
field and which are responsible for the scattering. 


C. One needs an understanding of the fundamental hydrodynamic mechanisms. 

In the current contribution we shall confine ourselves to problem area A. the inter- 
action between electromagnetic waves and an irregular ocean surface. (See also 
Gjessing, Hjelmstad, Lund (1984)). 


2 SCATTERING FROM A ROUGH SURFACE 
2.1 A summary of basic theory 

In recent contributions (Gjessing, 1981 a and b) , a somewhat detailed discussion was 
presented of the properties of the electromagnetic field scattered from a rough sur- 
face. 

The basic scattering equation can be expressed as 


E (K,t) ~ / f(r,t) e~ 2K ‘ T d 3 r (2.1) 

* —00 

-f — ► i 1 4 tt 

K = k . - k and |k| = — sin 9/2 

j S A 

A is wavelength and 0 scattering angle. 

In order to obtain backscatter from a scattering surface, the surface must have a 
structure, regular or irregular, containing scale sizes l = — ^ = A/2. 

Jv 

Referring again to equation (2.1), r is a position vector, t is time and f(r) is the 
delay function characterizing the surface. This complex function tells us how the 
scattering elements contributing to the "bulk scattering cross-section a" are distri- 
buted spatially. 

Resolving r and K into orthogonal coordinates in the horizontal plane 



E (K) - { / f(z) e Z dz} 6(K x x) (2.2) 

“00 

where 6 denotes a Dirac delta function. 

The scattered field B (K) vanishes unless K* is zero. Hence, E (K) * 0 unless & 
is normal to the wave crest. 

Let us assume that we illuminate the sea surface with a beam whose width is limited 
so as to cut out a section D of the ocean wave front (see figure 1) . 
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OCEAN WAVELENGTH. METERS 


Figure 1 


The geometry of the trunca- 
tion (convolution) process 


6 = 


-1 

tan 



Figure 2 The angular resolution of our 
radar varies with ocean wave- 
length. As the ocean wave- 
length decreases, the resolution 
approaches that of the antenna 
beamwidth. 


In this case E(K X ) will be of a sine form instead of the Dirac deltfA function 
when approximating the antenna beam function to a rectangular distribution. 

Thus, if D be the width of the truncation influencing f(x), and K = — be the 
wavenumber of the one-dimensional oeiay function f(z), (L is the ocean wa .^length) 
then the angular width (angle between nulls) 2a of the E(K X ) function is 



This is a relationship well known from antenna theory. The beamwidth 2a (azimuth 
resolution) is determined by the antenna aperture D expressed in wavelengths L. 
(See figure 1) . 

This is i lustrated in figure 2 and verified experimentally in figure 13. 

The radar antenna is illuminating a range interval B measured along the direction 
of wave propagation. This leads to a truncation of the delay function f(r). Then 

6F, « 0.16 £ for an exponential function (Gjessing 19Slb) . 

* B 

We shall now suggest n simple mathematical model baaed on basic and rather intuitive 
physical arguments to provide a simple way of assessing the experimental findings. 

We assume that long ocean waves (wavenumber approaching cut-off K Q ) are unidirection- 
al, The higher the wavenumber, the larger is the angular spread. For waves in the 
wavenumber region above K s , the irregularity structure is isotropic. We therefore 
write the complex irregularity spectrum simply as 

E(ic) = e(|k|,0) 

- K _n (cos 6 ) s 
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Figure 3 Complex wave spectrum cal- 


culated from the simple 
model presented 
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Figure 4 Scattering from turbulent 

irregularities caused by the 
"footprint" of a local wind 
field interacting with the 
wave structure of the sea 



where s, the planeness parameter, is expressed as 



where K Q is the wavenumber corresponding to the short wavenumber cut-off, and K g 
is the wavenumber at which the ocean surface structure becomes "random" and iso- 
tropic with no prevailing wave direction. 

A local wind field may through shear mechanisms break up the coherency of the ocean 
wave ridge. The effect of this is to widen the angular distribution of the back- 
scattered radar wave considerably# as illustrated in figure 4. 

Here we have obtained that the wave crests and troughs are broken up by turbulence 
and wave backing phenomena such that there is no "phase matching" between radio wave 
and ocean wave over the illuminated area contrary to the coherent case depicted in 
figure 1. 

Thus# when turning the antennas on azimuth angle 0# the radio wave sees a dominat- 
ing ocean irregularity scale which is increased from L when 0*0 to 

L * — - — 

0 cos 0 



This is illustrated in figure 4. Note that for certain values of K (and the cor- 
responding AF) the back-scattered radio signal will increase with increasing azimuth 
angle symmetrically on either side of the 0*0 direction. For experimental verifi- 
cation# the reader is referred to figures 11 and 13. 

It is well known from turbulence theory that the nearness of a rough boundary such 
as the sea surface results in a low wavenumber suppression of the turbulence field 
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and also to anisotropy (see e.g. H Panovsky 1960). (Modulation of the sea surface 
by the wind turbulence, anisotropic irregularity structure on the sea surface with 
an axis of symmetry determined by tha direction of the local mean wind (Gjessing 


1962) .) 


2.2 The multifrequency radar principle, frequency difference matching 

Up to this point we have Leen dealing with electromagnetic waves having wavelengths 
which are comparable with the irregularity scale of the scattering surface. 


We shall now investigate what information we can derive by matching some beat pa' 
terns (difference frequency) to the irregularity scales of the sea surface. Fig-’ 5 
illustrates the physics of the problem. 


-f 


±t 

6 GHz 


1 



For this AK matching scheme to work (bee 
e.g. Plant 1977), the scattering surface 
obviously must have irregularities which 
contribute to the delay function f(z) 
for scales corresponding to the difference 
wavenumber AK. If the gravity wave struc- 
ture shows up in the delay function f(z). 
i.e. if the spectrum of f(z) has Fourier 
components matching the wavenumber AK, we 
woul^ experience a radar return at the beat 
frequency Aw * cAK. 

There are many ostensible explanations for 
this "modulation" of the small scale struc- 
ture (to which the wavenumber K is matched) 
by the long wave phenomena (to which AK 
is matched). (See e.g. Plant 1977, PJ ant 
and Schuler 1980) . 


It is well known that the scattering cross- 
section of a surface increases with inci- 
dence angle. From a sinusoidal gravity 
wave we shall get maximum low angle return 
at the point of maximum slope. The turbu- 
lent intensity is governed by deunping factors such as shear forces (Lumley 1969) . 
Orbital motion within the wave structure changes the velocity shear periodically in 
phas^j with the gravity wave. Shadowing effects may enhance the patchiness of the 
radar illumination. Breaking waves and white capping resulting from the local wind 
fieJd etc complicates the issue still further. 


Figure c j The interaction of electro- 
magnetic waves with the sea 
surface 


Having presented intuitive arguments based on physical interpretations, now let us 
describe the principle of the multifrequency radar mathematically. 

Wfj start with the simple and general expression given earlier 

V(-) - E (K) - / f (z) d 3 z (2.4) 

C ' 

V 

-► m + •+ 

and we compute the covariance E (K) £ (K+AXj obtaining 

E*(K) E (K+AK) ' 9 ~3AK' z . I R (j) ."l** • d 3 r (2.5) 

R(r) is the complex autocorrelation of the delay function f(z). The phase factor 
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of the autocorrelation function _^(the term e Z ^) is rapidly oscillating with 

AK since the space coordinate z is large. 

_ -j 

The factor e represents the modulation of the spatial autocorrelation function 

by the radar carrier. Uut as the radar wavelength (5 cm) is orders of magnitude 
smaller than the scales of interest on the sea sur. face, this i./~tor will not affect 
the measured sea-surface signature because the individual scattei^rs move independ- 
ently such that the effect of the rapidly rotating phase factor will be averaged out. 

This can also be argued for mathematically by noting that the integral is of the 
form of a Fourier integral which is symmetrical with respect to K and AK. taking 
K and r as the Fourier pair of variables, and making use of the fact that the 
Fourier transform of a product is equal to the convolution of the Fourier transform 
of each factor, the integral converges to 

/R(J, . »!»)*♦<*•«> 

v . * (AK) (2.6) 

I 

wher <3 4 denotes convolution and 6 is Dirac's delta function. <X> ( AK) is the wave- 
number spectrum which again is the Fourier transform of the spatial autocorrelation 
function R(r). 

4 - -► -y 

As K N > AK, the term 0 (K) will vanish because K assumptionally is greater than 
the scattering medium cut-oif wavenumber. As an example, note that if £_ = 1000, 

(K) / <t> (AK) * 10 6 for a rectangular distribution of scattering elements.^ 


Normal izing this expression, we get the following familiar expression for a given 
fixed mcnostatic transmitter/receiver position 


E* (— ) E(^) 

~ c 


. A'jj . Au 

-3“ z , -3—' ' r 3 

e J R(r) e d r 


K>l ! 


(2.7) 


where A is a normalizing factor of the form / R(r) dr. 


Equation (2.7) states that the envelope of the complex correlation in the frequency 
domain of waves scattered back from a surface is given by the autocorrelation func- 
tion R(r) characterizing the surface. 


2.3 The motion pattern of the sea surface 

In the section above, we have shown that frequency components with mutual frequency 
spacing P couple to irregularity scale sizes AL * • 

We have already observed that the irregularity structure (the delay function) of the 
sea surface is determined by several mechanisms. Gravity waves are well behaved, 
and propagate at a velocity 


V - 




( 2 . 8 ) 


where L is the wavelength. 

Since the Doppler shift produced by a scatterer at velocity V is given by 

f m — jr.v 
1 2n 

- 2AF /^L 
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(2.9) 



giving 



We select a given AF^ (which is the same as selecting a given irregularity scale 
Lj), anJ we measure how the frequency covariance function 

R ( AF j , t ) * V(F,t) V*{ (F+AFj) ,t) 

is varying with time t. We then compute the power spectrum of the frequency covari- 
ance function R(AF^t). 

If then the ordered gravity waves dominate over the incoherent velocity components 
which ride on the gravity waves, we would expect the Doppler spectrum to have a 
maximum at the frequency corresponding to the d persion relation of gravity waves 
(eq (2,10) above), and we would expect a Doppler broadening determined by the velo- 
city spread 6V. 

3 EXPERIMENTAL RESULTS 

The adaptive multifrequency radar system developed by the authors 1 organization has 
been used for several applications: ship and aircraft identification, measurement 

of sea surface from a ground-based station and from an aircraft. A brief highlight- 
ing of results relaxed to the sea surface will now be given. 

3 . I Directional ocean wave spectra observed from the NASA Electra aircraft 

Illuminating the sea surface by a fixed 14° beamwidth side-looking antenna pointing 
10° downwards relative to the horizontal plane, the ocean wave spectra were deter- 
mined for various azimuth directions and for 15 different ocean wavelengths in the 
interval from some 16 o (couples to AF * 8 MHz) to 300 m (corresponding to AF * 

1/2 MHz). For each frequency separation AF the frequency covariance function 
V (F, t) V (F+AF,t) was computed. 

The NASA Electra operated north-east of Wallops Flight Facility (38.49.2 N - 74.24 W) 
at an altitude of 12 000 ft, January 20, 1983, time 2203Z. 

Banking the aircraft 10° so that the 4..itennas were pointing at a depression angle of 
20°, the aircraft completed a 360° circle. In this way the antenna was illuminating 
essentially the same area on the sea surface at all azimuth directions. The longitu- 
dina 1 dimension of the footprint is approximately 5600 m# whereas the transverse di- 
mension is approximately 2000 m. 

Flying in a closed circle# the frequency covariance function V(r) V*(F+AF) was 
computed for 15 different values of AF for all azimuth angles. The results are 
shown in figure 6. Note that the Bragg angle is clearly visible for the longer wave- 
lengths (215 - 88 m) . 

In terms of aircraft altitude H# depression angle e» and antenna beamwidth AB we 
get the following expression for the Bragg angle 

_ L sin a 

This is prsssntsd in figure 7. Note the good agrees. nt between theory end experiaen- 
tal results. 
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Figure 7 The measurement of Bragg angle 
is compared with theory. Note 
that the aircraf* altitude is 
12 000 ft, depression angle 30° 
and beamwidth 14 . 




Figure 6 



The aircraft flew in a closed 
circle pointing the antennas 
30° down from 12 000 ft. Jhe 
covariance function V(F)V (F+AF) 
is shown for 15 different 
values of AF. 
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Figure 8 Wave direction plotted a3 a 
fur ;ion of wavelength 
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Noting the two azimuth angles corres- 
ponding to destructive interference 
on either side of the wave direction, 
this can be determined with great 
precision as shown in figure 8. 

Finally, for an azimuth direction cor- 
responding to that of the propagation 
direction of the ocean waves, the dis- 
tribution with wavelength of the 
scattering cross-section is computed. 
This is shown in figure 9, Note that 
two wavelengths dominate, 130 m and 
18 m. 
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FREQUENCY SEPARATION (MHz) 
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Figure 9 There are two dominating 
wavelengths, 130 m and 
18 m. 
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OCEAN WAVELENGTH («) 

Figure 10 Spectra of wavelength and wave 
velocity. Note that wave- 
lengths shorter than that cor- 
responding to the spectral 
peak move at constant velocity. 
This is not the case for longer 
scales. The results are from 
Langesund 12 March 81 at 1030 
local time, the azimuth direc- 
tion being 13‘°. 


3.2 Directional spectra obtained from ground-based station on a cliff 

The complex field strength was recorded over a period of some *.2 minutes before the 
azimuth angle was changed, thus obtaining a plot of power spectral density and 
Doppler spectrum as a function of frequency separation (ocean wavelength) . A typi- 
cal result is shown in figure 10. 
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First, consider the results obtained on March 12, 1981 (1010 - 1235 local time). 

The weather and sea state can be characterized as follows: 

- Wind direction 90° (East) 

- Wind speed approximately 8 m/s 

- Wave height 1.4 m peak to peak 

- Waverider gives maximum wave height for f = 0.225 Hz 

- Wind decreasing at noon, turning northerly. Wave heignt reduces to 1 m peak to 
peak, choppy sea 

Figure 12 shows the wave intensity for a particular azimuth direction and also the 
Doppler spectra associated with various scale sizes. Note the very well defined 
spectral peak in the intensity versus wavelength distribution. Note also that the 
Doppler frequency is in very good agreement with the theoretical dispersion rela- 
tion. There seems to be a slight, but consistent, displacement to higher Doppler 
frequency. This is a result of an inward current. 
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Figure 11 Isolines of spectral intensity 
♦ (Ait) in the l^-K^ plane. 

Note the effect of changes in 
wind direction on the coherent 
irregularity spectrum. Note 
that the upper set of maps was 
obtained on 12 March 1981, 

1010 to 1135 local time, where- 
as the lower set refers to the 
situation during the time in- 
terval 1135 to 1235. 


Finally, in figure 11, we present the 
two-dimensional irregularity spectrum 
summing up che description of the sea- 
surface structure. Here we have plot- 
ted the isolines of spectral intensity, 
or morg specifically the quantity 
V (F) V*(F+AF)/V(F) 2 , in the K x ~Ky 
plane. This presentation should be 
compared with those of figures 2, 3 
and 4. Note also that the isolines 
of spectral intensity are ob- 

tained directly from a set of one- 
dimensional spectra such as those 
shown in figure 10. 

As will be seen, we dist- .guish bet- 
ween coherent and incoherent compo- 
nents. With "coherent component" we 
refer to the dominant Doppler spectral 
line caused by the "coherent" gravity 

wave (f = . 

The incoherent component is the "resi- 
due intensity" of the Doppler spectrum 
(the contribution from the skirts) 
obtained when the coherent contribu- 
tion wave is subtracted. It can be 
visualized as all the random turbulent 
velocity contributions outside the 
short velocity interval dominated by 
the gravity wave. Note, however, that 
although the velocity of the scatterers 
vary widely as seen e.g. from figure 10, 
the difference frequency AF selects a 
narrower range of irregularity scales 
centered round the scale size L given 
by 
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If these weak incoherent irregularities were entirely due to the local wind field, 
one would expect the local wind direction to provide an axis of symmetry. Unfortun- 
ately, the antenna systems were not flexible enough to be turned beyond the east 
direction (0 = 90°). It is therefore difficult to establish such a symmetry. The 
coherent contribution, however, has a pronounced axis of symmetry determined by the 
direction of the most energetic ocean wave. We see that the angular distribution is 
far from cos 2 0 distributed, long waves are more plane than the shorter ones, as sug- 
gested in figure 3. 


The lower set of directional spectra in figure 11 gives the situation one hour later 
when the wind has turned from east towards north. We see that the most energetic 
waves (the "eye" of the isoplot chart) have not changed direction, whereas the shor- 
ter scales and weaker disturbances are reoriented slightly. The distribution of the 
incoherent (turbulent) irregularities is drastically altered. 
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On 18 March 1981 there is no wind before 
noon. At 1330 the wind speed is 10 m/s 
from south, bringing the water height up 
to 2.5 m. The complex K plot shows a 
unidirectional spectrum. 

Two hours later, when the wave height has 
increased to 2.5 m peak to peak, and wind 
has turned easterly, the irregularity 
structure is very much altered. This most 
striking feature has (as indicated in 
figure 12 and also summarized in figure 13) 
the very pronounced appearance of a second 
peak at a wavelength of 25 m whereas the 
dominating spectral peak is still 50 m. 

Note that this two-peak spatial spectrum 
for a particular azimuth direction is not 
very pronounced in the wave-rider spectrum. 
This, presumably is a result of the omni- 
directional response pattern of the wave- 
rider. Note also that figure 12 shows a 
very pronounced two-peak structure in the 
Doppler spectrum for an ocean wavelength 
of some 70 m. Finally, figure 13 sums up 
the results <jf the afternoon run in the 
form of 0(AK) isolines for the coherent 
as well as for the incoherent components. 
Note that the wavelength corresponding to 
peak spectral intensity varies with azi- 
muth angle. 


Figure 12 Spatial and temporal wave 
spectra for an azimuth 
direction of 140°. Loca- 
tion is Langesundsf jorden, 
time is 18 March 81, 1525 
local time 
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Fj_gura 13 The two-dimensional ocean irregularity spectrum (isolines of $>(AK)) 

when the wav^ height has increased to 2.5 m peak to peak and the wind 
has turned towards east 
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